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We perform the quantum analysis of the light emitted by a synchronously pumped optical para- 
metric oscillator operating in the above threshold regime, i.e. when the peak power of pulsed 
pumping exceeds the threshold of continuous generation. We show that both regimes (below and 
above threshold) are realized at different times within each pulse of the signal field. We show that 
the quantum fluctuations of signal and pump pulses are correlated between nearby pulses at times 
that are placed in the same position relative to the center of the pulses, whereas fluctuations are 
totally not correlated at different times within the same pulse. The model also predicts the existence 
of cross-correlations between pump and signal pulses. 

It is shown theoretically that there is suppression of quantum noise at the frequencies multiple 
of the pulse repetition frequency in the spectra of phase quadratures of pump and signal fields 
measured by a balanced homodyne detection with a pulsed local oscillator. 

PACS numbers: 42.50.Dv, 42.50.Yj, 42.65. Re 
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Development and investigation of effective sources of non-classical multimode light are among current trends of 
quantum optics [H-Q ■ These sources appear as a main component required for parallel quantum information protocols 
0, @- From this point of view an optical parametric oscillator (0P0) pumped by a mode- locked laser with a period 
of pulses which is equal to the round-trip time of the oscillator cavity (so called synchronously pumped OPO or 
SPOPO Q) looks quite promising. A modal approach used in @, [1] showed that quantum state of signal light of a 
SPOPO operating below threshold in a degenerate regime generates a tensor product of squeezed states in different 
super-modes, each being a particular coherent superposition of longitudinal modes of different frequencies. Therefore 
the signal light generated by a SPOPO appears as an interesting resource for parallel quantum information protocols 
such as quantum teleportation and quantum key distribution. 

The same system was also considered in the time domain 0] using a simple physical model that neglects all 
dispersion effects in the parametric crystal. It was shown in particular that the quantum fluctuations of signal field 
are not correlated at different times within each individual pulse whereas they are correlated between nearby pulses 
at times that are placed in the same position relative to the center of the pulses. This paper dealt with the SPOPO 
operating below threshold. However it is well known that the light of a continuously pumped OPO also reveals 
nonclassical properties above the oscillation threshold. [Io| . Therefore the question about the quantum properties of 
SPOPO light above threshold looks quite natural. 

In the paper it was predicted that a multimode OPO (i.e. which is simultaneously resonant for several spatial or 
temporal modes) above threshold generates bright light in a single mode whereas the other modes stay in noncritically 
squeezed vacuum states, i.e. the amount of squeezing does not depend on threshold excess. The prediction was 
confirmed by an experiment on parametric generation of TEMio and TEMoi modes in squeezed vacuum states when 
the threshold is crossed. In our case it is interesting to check this prediction for the SPOPO, which is another 
multimode system. 

In this paper we will analyze the above threshold behaviour of SPOPO in the time domain where the system 
operation is described by two coupled equations for envelopes of pump and signal pulses. The physical model and the 
two-time technique used to describe the field evolution inside the oscillator cavity are similar to the ones applied in 
for the analysis of the below threshold regime. 

The paper is organized as follows. In Sec. [TT| we present our model of SPOPO and write Heisenberg-Langevin 
equations that describe in time domain its operation above the oscillation threshold. In Sec. IHII we estimate classical 
steady-state envelopes of generated pump and signal pulses. We assume that the pump pulses of oscillator have a 
rectangular shape. We determine quantum fluctuations of fields around the steady-state values in Sec. IIV1 For this 
purpose we solve linearized Heisenberg-Langevin equations applying adiabatic elimination of the pump field. As a 
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result we calculate and analyze pair correlation functions for the quadrature components of the output pulses. In 
the Sec. |VI] we consider the model of balanced homodyne detection of pulsed field and calculate the spectra of the 
photocurrent fluctuations. We generalize the results obtained for rectangular pump and local oscillator pulses to the 
case of arbitrary envelopes of pulses. 



II. PHYSICAL MODEL AND MAIN EQUATIONS 

A quantum analysis of the SPOPO below threshold on the basis of time approach has been performed in 0. As 
a result two-time correlation functions were obtained for the output signal field which correspond to the following 
picture of the system operation. Under pulsed laser pumping of a nonlinear crystal (see Fig. [T]), placed inside a 
ring optical cavity of the oscillator, pump photons are parametrically down-converted into pairs of correlated signal 
photons. The photons of each pair may leave the cavity in different pulses, during an overall time of the order of photon 
lifetime in a cavity giving rise to temporal correlations on the same range of time difference. Approximations 
of thin nonlinear crystal and instantaneous down-conversion, used in the work, lead to the peculiarity that field 
fluctuations at different times in a single pulse are not correlated. 

We will consider SPOPO operation above oscillation threshold on the basis of the physical model that we have used in 
the below threshold situation. Similarly we assume a degenerate parametric interaction for the carrier frequencies of the 
pump and signal modes, meaning that the following phase matching condition is fulfilled Afc = k p {uj p ) — 2k s {uj p / 2) = 0, 
where k p {u> p ) and k s (u) p /2) are wave vectors of carriers of pump and signal fields. The SPOPO ring cavity is resonant 
and of high-finesse both for pump and signal fields, so that the doubly-resonant configuration is realized. It takes the 
same time Tr for pump and signal pulses to make a single round trip inside the cavity and this time is equal to the 
pump repetition rate. We assume also that the cavity is dispersion-compensated by intracavity dispersive elements. 
This implies that optical pulses of arbitrary shapes are not distorted after one round trip inside the cavity. 

The above threshold regime is achieved by increasing the pump power. In this regime pump depletion becomes 
significant as well as the process of up-conversion of pairs of signal photons into pump photons. Under these conditions 
the system becomes nonlinear. In particular one can expect that correlations between the pump and signal pulses 
will develop. Therefore an adequate system description must consider the coupled evolutions of both fields. 

In order to describe system operation we use the time-domain approach presented in detail in Appendix [21 The 
Heisenberg-Langevin equations obtained under this approach have the following form [9( 



dA p (t, T) 
dT 



= -Kp (A p (t,T) - A (t)) -gA s 2 (t,T)+F p (t,T), (1) 
dA f T T) = -KsA a (t, T) + 2gA p (t, T)A\ (t, T) + F s (t, T) (2) 

Here A p {t,T) and A s (t,T) are quantized envelopes of pump and signal pulses inside the cavity. Time argument 
t is treated as the time deviation from the pulse center and it changes in the interval from — Tr/2 to Tr/2; the 
dependence of envelopes on the second time argument T describes their changes from pulse to pulse. The amplitudes 
are normalized so that values (A p s A PyS ) have the meaning of mean fluxes in photons per second through the cross 
sections of the light beams. k p and k s are the loss rates of pump and signal fields respectively; g is a constant 
characterizing the parametric coupling; A (t) is the classical steady-state envelope of the pump pulses inside the 
cavity, depending only on time t since the pump pulses are supposed to be perfectly identical. One can write this 
envelope in the following general form 



A (t) = 0VoW e*" (i) , (3) 

where Ng(t) is the intensity shape of the pump pulses and 4>i n {t) their possible phase modulation. 

The last terms in equations F p and F s describe Langevin noise sources which are related to the vacuum fluctuations 
of the incoming fields. In this case the fields have zero mean values and characterized by the following nonzero pair 
correlation functions (for details see Appendix [B]) 

(F r (t,T)F}(t',T')) =2K r 5(T-T')5(t-t'), r = s,p (4) 



III. OSCILLATION THRESHOLD AND SEMICLASSICAL STEADY-STATE SOLUTIONS 

As it is known (see, for example, [13]), & n optical parametric oscillator is a system that exhibits properties of a 
second-order phase transition when pump power achieves a certain threshold value. In our notation the threshold 
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power reads 



N th = (5) 

For pulsed pumping the threshold is achieved when the peak power of pulses is equal to this value. Evidently, 
corresponding mean power can be much less than the threshold value. This is an important advantage of pulsed 
pumping with respect to the continuous one. 

It is important to note that in the Heisenberg-Langevin equations ([1} and {2J time t appears as a parameter and 
the Langevin noise sources are delta-correlated relative to this time according to (U]). Hence in our model different 
temporal parts of an individual pulse (pump or signal one) develop in time T independently of other parts with 
corresponding instantaneous value of the external pump power No{t). As a result both regimes (below and above 
threshold) could be realized within the single pulse at different time instants if peak pump power exceeds threshold 
value (see Fig. [2]). This situation is analogous to what happens in the spatial domain in a degenerate confocal 
cavity [l^. Since the properties of pulsed fields below threshold were analyzed in Q we will consider further the 
simplest case, when external pump pulses have rectangular temporal shape of duration t p with intensity exceeding 
threshold one. In the last section we will generalize our results to arbitrary envelopes of the pulses. 

Let us define a time dependent pump parameter: 



For rectangular pump pulses the parameter is non-zero no > 1 when t changes from — t p /2 to t p /2. Between pump 
pulses the generated fields are vacuum noise while within the pulses they develop in a way that is described by coupled 
Heisenberg-Langevin equations. At this time intervals the fields are characterized by non zero average amplitudes. 
Here we determine the average steady-state amplitudes denoted as A p and A s . In the next section we will analyze 
the properties of the fluctuations around the average The amplitudes must satisfy equations which are the classical 
counterparts of Heisenberg-Langevin equations ([TJ and @. Technically, we replace the operator quantities in these 
equations by the c-number quantities and to reject the noise sources. Finally setting partial derivatives equal to zero, 
we find the stationary solutions in the form 

A p (t) = yfT p e ^"W, A s {t) = ±^fN~ s ( 7 ) 

Classical intensities N p and N s are defined by the following expressions 

N p = N t h , N s = ^ (mo - l)N th , (8) 

K s 

which agree with the fact that the average amplitude of pumping field inside the cavity (and respectively at the 
output) does not depend on the pump intensity [13j. Also one can see that there are two possible opposite phases of 
the signal field for a fixed pump phase. One can show that both solutions are stable jioj . 



IV. SOLUTIONS OF LINEARIZED HEISENBERG-LANGEVIN EQUATIONS 



The following expressions define the quantum fluctuations of envelopes of pump SA p and signal 5A S fields within 
the pulses (i.e. at time interval — t p /2 < t < t p /2) 



(9) 
(10) 



A s (t,T) = {±^N~ s + 5A s (t,T))e 1 ^' 2 

Let us assume that the fluctuations are small compared to the mean values 

5A r (t,T) < y^Vr. r=p,s (11) 

It will be shown later that above threshold these inequalities are well-satisfied. Substituting expressions © and (TTU)) 
into (fT]) and ([5} and ignoring the second-order terms according to the assumption of small fluctuations, one obtains 
the following linearized equations for the fluctuations of fields 

d ; 5A p (t, T) = -K p 6A p (t, T) T- 2g^N~ s 5A s {t, T) + F p {t, T), (12) 



dT 
d_ 
df 



SA s (t, T) = -K s SA s (t, T) + 2g^W p 5A\(t, T) ± 2g^N~ s SA p (t, T) + F s (t, T) 



(13) 
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These equations can be conveniently solved by writing them in terms of quadrature components of fields defined as 
real and imaginary parts of complex amplitudes: 



SA r (t,T) = 5X r (t,T) + i6Y r (t,T), 



(14) 



where SX r — 5X]. and SY r — SY^ . This choice of quadrature components is such that their fluctuations define 
fluctuations of intensities and phases of fields, respectively. These quadrature components can be experimentally 
measured using balanced homodync detection of fields under appropriate choice of the local oscillator field that will 
be considered in the Sec. I VII 

One gets the following equations for quadrature components of pumping field 

° ■ ~" -K p 5X p (t,T) T 2g^N s SX s (t,T) + F p (t,T), 



dT 
d_ 

dT 



5X p {t,T) 



and signal field 



d_ 

dT 
0_ 

dT 



SY p (t,T) = - Kp 6Y p {t,T)T2gVN~sSY s {t,T)+F£(t,T), 

6X s (t,T)=±2gy/N~ s 5X p (t,T) + F&,T), 
5Y s (t, T) = -2n s SY s (t, T) ± 2g^/N~ 5Y p (t, T) + F'J{t, T) 



(15) 
(16) 

(17) 
(18) 



Here we took into account that above threshold N p — K^/(Ag 2 ) and we also defined hermitian quadrature components 
of Langevin sources: 

F p (t, T)e- 4 *«" = F p (t, T) + iP;'(t, T), 
F s (t,T)e-^" /2 = F' a (t,T) + iF'J{t,T), 
with the following correlation functions, obtained from (jU): 



(19) 
(20) 



&{t,T)FW ,T')> = (F?(t,T)F?(t',T')) 



S(T-T')S(t-t'), 



(F;(t,T)F;'(t',T'))^-(F;'(t,T)F;(t',T'))= i ^-S(T-T')5(t-t / ) 



(21) 
(22) 



Expression (|22p shows that the quadratures of fields are correlated. Thus fluctuations of X- and Y-quadratures of 
pump and signal fields will be also correlated. Nevertheless in the case we are interested in these correlations are not 
important. 

Let us consider the case when relaxation of the pump field in the cavity with the rate k v is the fastest process. 
Then equations (|15[) -(|18 p can be solved using adiabatic elimination of this field. Setting derivatives equal to zero in 
equations (fT5]) - (fT6|) . the fluctuations of the pump field have the form 



5X p (t,T) 
8Y p (t,T) 



/ 

Ta/ — 

Kp 



5X s (t,T) 
SY.(t,T) 



1 



F"(t,T) 



(23) 



Inserting these expressions into equations (fT7|) and (fl8|) . one gets simple differential equations which have the following 
solutions 



6X s (t,T) 
SY s (t,T) 



dT' 



± [«*( F p^ T ') 
V«p \F p '(t,T') 



F'J(t,T) 



,{T-T') 



(24) 



where k x — 2k s (/xo — 1) and n y — 2k s /xo are effective damping rates of fluctuations. According to adiabatic approx- 
imation these rates fulfill the condition n x , k v <C k p that limits pump parameter /io k p /k s . We will show further 
that the main quantum effects appear close to the oscillation threshold when fj,o ~ 1, which is in agreement with our 
restriction. 

These solutions enable us to determine the properties of the output fields which are of practical interest. For this 
purpose one uses the boundary condition on the output mirror that reads 



A° ut (t,T) = y/%Ar(t,T) - y/T=%4?(t,T). 



(25) 



Here T T is the transmission coefficient of the cavity mirror related to the loss rate of the field K r = T r /(2Tn) (when 
%• <C 1). Also vacuum fluctuations of the incoming field A 1 ™ in the expression are related to Langevin noise sources 
F r according to (|B1|) . The boundary condition is also valid for the quadrature components of fields. 
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V. CORRELATIONS BETWEEN PULSES 

Let us now determine the correlation functions of the quadratures of output fields. Before it we turn from time 
t that describes a deviation from center of pulses to usual time scale. For this one returns to the discrete pulse 
numbering, replacing time T with discrete number n and time t with t — nT R in the following way: 



T — » nT R , t^t-nT R , T R S(T - T') -> 6 nn > , 

SX™\t, T) -> SXZ\t - nT R ), 5Y° ut (t, T) -> 5Y°^{t - nT R ). (26) 



Thus, using results of the previous section, we obtain the following pair correlation for the quadrature components of 
pump pulses 

(SXfflt - nT R )5X° p %{t' - n'T R )) = ± (6 nn ,S(t - t') + 2 Ks T R e- 2K ° T ^°~^ n - n '^5(t - t' - (n - n')T R )) , (27) 
(ri£f(t - nT R )6Y£$tf - n'T R )) = ~ (s nn ,5(t - t') - 2k s T r 1} e -^sT R ^\n-n'\ 6{t _ t > _ (n _ ^ )Tr) ^ (2g) 

and signal pulses 

rt£(t - nT R )5X°$(t' - n'T R )) = ± ^ w <J(i - t') + -^- e -2^ R ( w -i)|„-„'| 5(i _ ^ _ (n _ „/ )Tfl ^ ^ (29) 

(«5Y™'(i - nT R )5??%V - n'Tfl)) = 1 (<W<^ - t') - ^A e -^ R ^\ n -W\ s ^ _ t ' _ ( n _ n ' )Tfl) ') . (30 ) 

4 V Mo / 

These expressions are similar to the correlation functions of pulsed signal field generated by SPOPO below oscillation 
threshold. The first terms on the right hand side which are proportional to 5(t — t') are due to the incoming vacuum 
field reflected by the coupling mirror of the cavity. The second terms are related to the fields coming out of the cavity. 
For individual pulses, i.e. when n = n', one can neglect second terms which are proportional to k s T r <C 1. Thus the 
fluctuations of the quadratures in the pulses are vacuum fluctuations, except the fluctuations of X-quadrature which 
become infinitely large close to threshold when —> I. Since these solutions are obtained under the assumption of 
small fluctuations one must define how close can we approach the threshold from above. It is shown in the Appendix [C] 
that for typical experimental parameters the condition of small fluctuations is fulfilled and the solutions are correct if 
(/xo-l) > 10~ 3 . Therefore we will use the following minimal value of the pump parameter in quantitative estimations 
Mo = 1-1. 

There are also quantum correlations between different pulses (when n ^ n') for the X- quadrature of the fields, and 
anticorrelations for the Y-quadrature with the correlation coefficient proportional to k s T r . The number of significantly 
correlated successive pulses can be evaluated by the factor in the exponential, which is proportional to (k s T r ) . The 
delta- function shows that the correlations between different pulses have a "local" character: they are effective only 
when the time differences are a multiple of the period T R . 

Correlation coefficient and number of correlated pulses also depend on pump parameter fiQ (i.e. on pump power). 
For Y-quadrature of the signal field they are maximal close to threshold, when /xq ~ 1- For Y-quadrature of the 
pumping field correlation coefficient tends to zero in the vicinity of threshold. Away from threshold the correlation 
increases while the number of correlated pulses decreases. Also close to threshold the fluctuations of the X-quadrature 
of signal field becomes infinitely large, which restricts the minimal value of pump parameter as remarked above. 

Solutions (|2"3")l and (|24[) also show that fluctuations of pump and signal fields are correlated in contrast to below 
threshold regime. One can get the following symmetrized cross-correlation functions for these fields at the output of 
the oscillator 

{5X° p %{t - nT R )8X° s %{t' - n'T R )) + (8X%(t - nT R )5X° p %,{t' - n'T R )) = 

" KsTr ^ 2{^- 1) e- 2K ' T ^-^ n - n '\8{t -t'- (n - n')T R ), (31) 

(6Y p °X(t - nT R )5Y$tf - n'T R )) + (5Y^{t - nT R )5Y p °$(t' - n'T R )) = 

~ KsTr J it^l e ~^T RfM An-n'\ S{t _ ^ _ („ _ u ' )Tr} (32) 
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The properties of the cross-correlations (number of correlated pulses, correlation coefficient and local character of 
correlations) are analogous to the properties of fields itself that we have considered above. Close to the oscillation 
threshold the correlations between pump and signal fields become infinitely large for X-quadrature and decrease for 
Y-quadrature. The minus sign indicates that the fluctuations are anti-correlated. 



VI. QUANTUM EFFECTS IN THE SPECTRA OF FIELDS 



Let us remind that the temporal features associated with the correlations of pulses are small (of the order of 
k s T r <C 1), while the number of correlated pulses is defined by inverse value (k s T r ) , Therefore one expects 
to obtain important quantum effects by observing integral characteristic of fields, such as the Fourier spectrum of 
fluctuations of quadratures. Let us consider the measurement of field quadratures obtained by a balanced homodyne 
detection of the output fields (see Fig. [3j. The signal or pump field is mixed at a symmetric beamsplitter with an 
intense local oscillator field of the same optical frequency (u> p /2 or u) p ). In this case the fluctuations of the difference 
photocurrent are given by the expression 

Si r (t) = /3*(t) 5A° ut (t) + 0(t) 5A° ut Ht), (33) 

where f3(t) is the complex amplitude of the local oscillator. We suppose here that the local oscillator field is a train 
of identical pulses and that their period is equal to the period of analyzed pulses T R 

p(jb) = Y,Po(t-nT R ). (34) 

n 

The envelope of pulses has the form 



ftW = V^e i W' + *) 1 (35) 

where phase modulation of pulses is matched with the detected field: for pump field tp(t) — <j>i n (t), for signal field 
ip(t) = <f>i n (t)/2. As a result choosing constant phase shift $ = or $ = 7r/2 one finds: 

Si r (t) = 2 VN L o(t n ) , t n = t- nT R . (36) 

The shape of local oscillator pulses, their duration and delay relative to the pulses of the fields to analyze can be 
experimentally adjusted. We suppose here that the pulses have a rectangular shape of duration r p and are ideally 
synchronized with the pulses of analyzed field. In this case the current fluctuations are equal to the quantum 
fluctuations of the fields inside pulses. In Sec. IVIII we will generalize our results to arbitrary envelopes of pump and 
local oscillator pulses. 

Let us now consider the measurement of Y-quadratures of pump and signal fields. Substituting expressions 
and pop in (1361) . we derive the pair correlation functions for the currents: 



(5i(t) 5i(t')) p = V N LO (t - nT R ) [5(t - t') - 2k s T r ^° H e - 2 ^o|t-t'| § (< _ t > _ ( 
(5i(t) Si(t')) s = N Lo{t - nT R ) U(t - - £ 5{jb _ f _ (n _ n ')T R )\ 

n \ ^ n' ) 



n - n')T R ) , (37) 



(38) 



One sees that the correlations between pulses lead to temporal periodic correlations of the photocurrent. Let us 
determine the frequency spectrum of the photocurrent quantum noise defined as: 

+T/2 +T/2 

(& 2 ) n = r lim ^ j dt J dt' (Si(t) 6i(t')) e m( f ' (39) 

-T/2 -T/2 

Substituting (j3"T|) and ([3"5]) into (|3T)1) one gets the following explicit expressions for the spectrum of the detected pump 
field 

(Kh =l _ 8^(^0-1) (40) 

(I) m it 2 ... + (fi _ 2nm/T R f ' 



7 



and signal field 

(I) m= tt, tetfi + (Sl- 2nm/T R ) 



Here fluctuations are normalized to the corresponding shot (quantum) noise of local oscillator 

+Tr/2 

W = jr J dtN Lo(t), (42) 

-Tr/2 

which for rectangular pulses is equal to (/) = (t p /Tr)Nlo- 

These expressions show that the shot noise is reduced around all resonant frequencies of the cavity VL m — 2itm/T R . 
This coincides with the result obtained in the below threshold regime. For pump field a maximal noise reduction of 
fifty percent is achieved at /io — 2. For signal field the expression predicts full reduction of the noise at the oscillation 
threshold, i.e. (Sil)n m /(I) — > when /x — > 1. 



VII. ARBITRARY SHAPES OF PULSES 



Our results can be easily generalized to arbitrary shapes of pump pulses with corresponding pump parameter 
fi(t). This is possible because in our model different temporal parts of individual pulses develop independently of 
each other inside SPOPO. Hence at time instants when instantaneous pump power exceeds parametric threshold, 
i.e. fj,(t) > 1, the average amplitude and fluctuations of signal field are defined by d)-© and ([2^|) - (|5tl|) . respectively, 
where replacement /iq — > fi(t) is made. The average amplitude of intracavity pump field is stabilized at the threshold 
value according to G])-© and its fluctuations are characterized by correlations (|2"T|) and (|2"5|) . At other time instants, 
when /j,(t) < 1, the signal field is in the below threshold regime. 

Also it is not difficult to generalize the results of balanced homodyne detection, namely the spectra (|4"0)l and (gT| , 
that we have derived for rectangular pump and local oscillator pulses. Let us assume that the local oscillator pulses 
with the envelope Ni,o{t) are quite short (tlo < T p) and synchronized with the analyzed pulses to probe temporal 
parts of the pulses generated above threshold. Then general expressions can be obtained by averaging intensity shape 
of the local oscillator pulses with expressions and (|4ip . where pump parameter depends on time fi(t) 



+Tr/2 



(Kh = Y- f dtN LO (t)(l- £ 

R J ,„ V 771=0,1,5 

+Tr/2 . 

(Sil)a = Y- f dtN LO (t)ll- 

R J ,„ \ 771=0,1,5 



lK 2 Mt) - 1) 



2 



>2 



(43) 



(44) 



These expressions mean that current fluctuations at a particular frequency appear as a weighted sum (integral) 
of fluctuations from all non-correlated parts of the measured pulses. Consequently, the detected quantum noise is 
sensitive to the temporal properties of the local oscillator pulses, particularly to their duration and to the delay 
relative to analyzed pulses, as in the below threshold regime. Fig. [4] presents suppression of the photocurrent shot 
noise at zero frequency as a function of delay At of short local oscillator pulses {tlo "C t p ) relative to the signal ones. 
Curves are obtained for Gaussian pump pulses of SPOPO fi(t) = ^e^ 2 ^^^ 2 and with different values of peak power, 
including below threshold values. In order to determine the noise suppression at the edges of the signal pulses, where 
the field is in the below threshold regime, we used the following expression of the photocurrent spectrum from Q 




* / dt N LO (t) 1- Y _]. (45) 



a T „ \ ^ ,2 



This expression was also used in the case when the peak pump power is less than the threshold value, i.e. fiQ < 1. 

The figure shows that the fluctuations of signal field are nonuniform inside pulses. In the central part of bright 
pulses (/io > 1), where the field intensity is maximal, the fluctuations are larger than at the edges. Of course measuring 
local squeezing with the help of infinitely short LO pulses is a purely theoretical scheme. In the realistic case when 
pulses of local oscillator have finite duration, the photocurrent fluctuations must be weighted by the intensity of the 
local oscillator pulse that corresponds to expression (|4~4"1) . 
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VIII. CONCLUSION 



Let us summarize our results: if the peak pump power exceeds the threshold of continuous OPO oscillation then 
both regimes (below and above threshold) coexist at different times within the pulses of signal field. Bright parts of 
the pulses above threshold are characterized by an average amplitude. At the edges of pulses the field has properties 
of the below threshold case. The quantum fluctuations of pulses turn out to be totally not correlated at different 
times within the same pulse, whereas they are correlated between nearby pulses at times that are placed in the same 
position relative to the center of the pulses. Above threshold the model also predicts the existence of correlations 
between pump pulses with the same properties and cross-correlations between pump and signal pulses. 

These correlations can be measured in the scheme of balanced homodyne detection using a pulsed local oscillator 
synchronized with the required times of the analyzed field. We have also shown that these correlations lead to 
a suppression of quantum noise in the spectra of phase quadratures of pump and signal fields around frequencies 
Q m = 2Trm/Tn (m = 0, 1, 2, . . .), where Tr period of pulses. A stronger noise suppression is achieved under detection 
of signal field when thelocal oscillator pulses are delayed relative to the peak of signal ones. 

Since a SPOPO is a multimode system it is interesting to compare our results with the prediction of the paper pl| . 
According to this paper, when one increases the pump power, the multimode parametric oscillator starts oscillating 
in the mode with the lowest oscillation threshold whereas other modes stay in non-critically squeezed vacuum states. 

One can consider for SPOPO that the bright parts of the signal pulses correspond to the modes developing above 
threshold and the field at the edges of pulses is formed by the below threshold modes. The correlation properties of 
the field show that an individual mode is a train of correlated delta pulses whereas different modes appear as trains 
of pulses delayed to each other. 

It is well-known that the balanced homodyne detection technique detects optical field amplitudes in a particular 
spatial-temporal mode defined by a coherent local oscillator pulse [lj| . Therefore shot noise suppression at the Fig. 0] 
could be treated as a squeezing of the stated modes of signal field since different delays of the local oscillator correspond 
to the detection of different modes. Further interpretation of the results in terms of the modes apparently should take 
into account dispersion effects in the parametric crystal such as mismatch and dispersion of group velocities. 

The study was performed within the framework of the Russian- French Cooperation Program " Lasers and Advanced 
Optical Information Technologies" and the European Project HIDEAS (grant No. 221906). It was also supported by 
RFB-R (No. 08-02-92504). 
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Appendix A 

In the degenerate parametric generation configuration that we consider here the field operator inside the cavity is 
equal to 

E(z,t)=E p (z,t) + E s (z,t). (Al) 

We use the plane wave approximation, so that the field amplitudes depend only on one longitudinal coordinate z 
measured along optical axis of the cavity. As the pump of the parametric crystal is realized by a train of optical 
pulses of duration close to 100/s it is possible to disjoint in the standard form quick oscillations of fields with optical 
frequencies uj PiS from slow changes of their envelopes [l5j . At the crystal entrance (z = 0) for the pump p and signal 
s waves the two fields read 

E r (0,t)=i( — 2 cS ) e~ lWrt M0,t), r = s,p, (uj p = 2u„), (A2) 

where n r = n r (uj r ) are the indices of refraction of the crystal. Following expression takes into account the periodic 
temporal structure of the fields 

A r {0,t) =J2 A r,n{t-nT R ). (A3) 

a 

Here A r ^ n (t — nT R ) is the envelope of the n-th pulse. Argument of the envelope t — nT R describes time deviation from 
the pulse center and changes in the interval from —T R /2 to T R /2. 

In order to describe evolution of fields inside the SPOPO cavity we use the two-time approach applied by Haus 
in Ref. [l6[ for developing a quantum theory of actively mode-locked lasers. We assume that the envelopes of pulses 
are not significantly changed from one pulse to the next, an hypothesis that is typically valid in experiments with 
high-finesse cavity and weak parametric amplification. Then the dependence on discrete number n could be replaced 
approximately by the a continuously varying temporal parameter T in the following way 

A r!n (t-nT R )^A r {t,T), (A4) 

where t on the right hand side of the expression denotes time deviation from the pulse center. Thus expressions (|A3|) 
and (|A4[) esteblish correspondence between envelope of the field A r (0,t) and envelope of an individual pulse A r (t,T) 
that depends on two time parameters. 



Appendix B 

Let us determine properties of Langevine noise sources entering in the Heisenberg-Langevine equations HJ and © . 
The sources are caused by vacuum fluctuations of incoming coherent /vacuum fields. Following expression directly 
connects these quantities 



F r (t) = ^-A™%t) « ^A" r ac (t), r = s,p, (Bl) 

where T r is the transmission coefficient of the mirror at the corresponding frequency. Vacuum fluctuations are 
characterized by the following non-zero correlation functions 

(A v r ac (t)A v r ac t(i')) = S(t - t'). (B2) 

Using written expressions and pulse representation of slow envelopes (|A3|) one gets that noise terms should satisfy 
following relations 

(F r , n (t - nT R )F} n ,(t' - n'T R )) = ^5 nn ,S(t - nT R - {t 1 - n'T R )). (B3) 

Finally making transition to continuous time parameter (|A4p and using relations S nn > — > T R S(T — T') and t — nT R — > t, 
one gets required correlators 
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Appendix C 



Here we determine values of pump parameter at which the condition of small fluctuations is fulfilled: 6A r (t,T) <C 
\/N r . It was shown in Sec. [V] that above threshold fluctuations of fields are vacuum except X-quadrature of signal 
field whose fluctuations become infinitely large in the vicinity of threshold when /j,q 1. Therefore condition of small 
fluctuations could be written in the form 

(6X*) « N s , (CI) 

where (SXg) is a variance of X-quadrature of intracavity signal field. In order to determine it let us evaluate correlation 
function of the quadrature (SX s (t, T) 8X s (t' ', T')). The function could be easily obtained using explicit expression for 
fluctuations of X-quadrature and correlation functions for Langevine noise sources (|2ip . One gets 

(6X s (t, T) SX s (t', 7")) = \ (l + 2(^-1) ) e- 2 "-<"°- 1 >< r - r '>a(t - if) (C2) 

We have to put T — T' and t = t' in the expression in order to find quadrature variance. However divergence 
appears due to the fact that field fluctuations are not correlated in each individual pulse. One can avoid the difficulty 
defining dispersion in the following way 

t+T F /2 

{5X 2 S ) = ±- f dt' (SX s (t,T)SX s (t',T)) (C3) 



t-Tp/2 

where Tp is an averaging time. Thus dispersion is given by 

K i+ ^) (C4) 

Finally let us take into account that average classical intensity of signal field tends to zero in the vicinity of threshold 
according to expression N s = (po — l)N t h- Therefore the condition of small fluctuations is fulfilled if: 




In order to estimate this restriction let us choose following parameters Q: pump wavelength A p = 0.4/im, threshold 
power of continuous generation P t u — 50W, averaging time (which is equal to the correlation time) Tp = lOfs and 
k p = 10 k s . As a result one gets 

Mo - 1 » 10~ 3 . (C6) 
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FIG. 1: Synchronously pumped optical parametric oscillator (SPOPO): (a) pulsed pumping field; (b) parametric down- 
conversion of pump photon in a pair of signal photons inside a nonlinear crystal; (c) establishing of quantum correlations 
between pulses of output signal field. Solid line - pump field; dashed line - signal field. 




(b) 

FIG. 2: Two regimes of oscillation within an individual pump pulse of SPOPO: (a) above threshold regime; (b) below threshold 
regime. Nth - threshold power of continuous generation; No(t) - instantaneous power of pulsed pumping. 
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FIG. 4: Shot noise reduction of difference photocurrent at zero frequency in dependence of delay At of short local oscillator 
pulses (tlo <S t p ) relative to signal ones. Pump pulses are Gaussian with different values of peak power /xo- 



